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Many physical systems respond to slowly changing external force through avalanches spanning broad range 
of sizes. Some systems crackle even without apparent external force, such as bursts of neuronal activity or 
charge transfer avalanches in 2D molecular layers. Advanced development of theoretical models describing 
disorder-induced critical phenomena calls for experiments probing the dynamics upon tuneable disorder. 
Here we show that isomeric structural transitions in 2D organic self- assembled monolayer (SAM) exhibit 
critical dynamics with experimentally tuneable disorder. The system consists of field effect transistor 
coupled through SAM to illuminated semiconducting nanocrystals (NCs). Charges photoinduced in NCs 
are transferred through SAM to the transistor surface and modulate its conductivity. Avalanches of isomeric 
structural transitions are revealed by measuring the current noise 1(f) of the transistor. Accumulated surface 
traps charges reduce dipole moments of the molecules, decrease their coupling, and thus decrease the critical 
disorder of the SAM enabling its tuning during experiments. 

Many physical systems respond to slowly changing external force through discrete, impulsive avalanches 
spanning a broad range of sizes 1 " 3 . A common feature in such systems is that statistics of various 
avalanche parameters, height, size, or duration, lack a characteristic scale. The far from complete list 
of systems exhibiting crackling avalanche dynamics includes earthquakes and landslides, Barkhausen noise in 
ferromagnets, crumpling of paper wraps, acoustic emission in martensitic transformations, fractures propaga- 
tion, popping foam bubbles, flux penetration into superconductors, and many others 1 " 3 . There are also systems 
that crackle without application of apparent external force, such as bursts of neuronal activity in cortex tissue 4 or 
charge transfer avalanches in 2D molecular layers 5 . Advanced development of theoretical models describing 
disorder-induced critical phenomena calls for experiments probing the dynamics upon tuneable disorder 1 " 3 ' 6 . 
Scale invariance in systems with avalanche dynamics is often attributed to the proximity of a non-equilibrium 



phase transition 1 3 . Avalanche sizes, S = 



I(t)dt and durations, T, in the vicinity of critical disorder R c have power 



law distributions D(S) oc S -T and D(T)ccT~ a , with universal critical exponents t and a depending on the dimen- 
sionality and type of criticality of the system. These exponents are related by the critical exponent relation 
a — 1 = (t — 1)/ crvz 1 " 3 . The product avz determines the exponent of the noise spectral density, which is expected 
to scale at high frequencies as S(co)ccco~ 1 ^ <TVZ and appears also as the power law exponent that characterizes the 
dependence of the average avalanche size <S> on the avalanche duration (S) oc T l ' avz . 

Along with exponent relations the most characteristic feature of the critical dynamics is the existence of various 
scaling collapses involving detail-independent (universal) critical exponents and scaling functions 1 " 3 . 

The detail dependent (non-universal) critical disorder R c is always determined in terms of the coupling 
strength in the system /. For the 3D random field Ising model RJJ = 2.16. For infinite 2D systems R c = 0, 
irrespective of the coupling strength. However, space limited 2D systems scale as if critical disorder is greater than 
zero and the apparent ,R C becomes system size dependent 6 . Disorder in the popular random field Ising model is 
represented by the width of distribution of quenched random magnetic field 1 " 3 . In strongly ordered systems, R <C 
R c , the local field is fairly uniform and a single spin flip can trigger the entire system to flip. With increasing 
disorder, at R = R c the spin flip avalanches become scale invariant and cover the entire range of sizes. With disorder 
increasing slightly above R c the avalanche size distributions show a disorder dependent cut-off 1 " 3,6 . However, 
experimental verifications of effects of changing disorder on the dynamic properties of avalanches are still missing. 

In this paper we report on a simple molecular NCs system in which the critical disorder can be in- situ tuned by 
light illumination. We demonstrate that such disorder changes can be probed by monitoring charge transfer 
dynamics through the monolayer. 
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Figure 1 | Schematic drawing of the experimental device. 

Charge transfer through organic molecular frameworks was 
extensively studied both theoretically and experimentally 7 " 9 . Quan- 
tum tunnelling dominates conductivity of thiolated alkenes 10 . The 
tunnelling barrier depends on the type and length of molecules and 
on their isomeric structure 11 . The barrier for all-trans state is much 
lower than the barrier for the gauche state 12 . We have recently shown 
that charge transfer within a molecular layer may have a collective 
avalanche character associated with disorder driven criticality and 
nonequilibrium phase transition 5 . 

Results 

Schematics of our experimental device are shown in Fig. 1. Colloidal 
InAs/CdSe/ZnSe nanocrystals (NCs) are covalently absorbed on a 
self- assembled organic monolayer that links them to the surface of a 
two dimensional electron gas (2DEG) GaAs/AlGaAs field effect tran- 
sistor (FET). SAM linker of 1,9 nonanedithiol HS-(CH 2 ) 9 -SH 
(NDT), or 1,4 butanedithiol HS-(CH 2 ) 4 -SH (BDT) molecules is 
self-deposited from an absolute ethanol solution 13 . The device is 
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Figure 2 | Avalanche duration histogram (a), size histogram (b), and the 
average avalanche size as a function of the duration (c), calculated from 
the first measurement for NDT device under illumination with a 0.2 dB 
attenuator. Solid lines represent best fits to power low distributions. 



illuminated by light with energies smaller than the GaAs energy band 
gap and larger than the InAs energy gap. When the NCs absorb light, 
the photo -excited holes are transferred to the FET surface states via 
SAM, change the surface potential, bend the conduction band, 
increase the 2DEG electron density, and thus increase the conduc- 
tivity of the transistor channel 14 . This has been confirmed by mea- 
surements of the DC response of our system to light illumination, 
showing monotonous decrease of the channel resistivity with illu- 
mination time. A change in the FET current can be orders of mag- 
nitude larger than the number of photo electrons excited in the 
NCs 15 . 

Fluctuations in the FET conductivity reflect directly the dynamics 
of charge transfer through the SAM layer. When collective charge 
transfer appears, the FET noise is dominated by large non-Gaussian 
avalanche-like fluctuations. The trigger that initiates an avalanche is 
a spontaneous transition of a molecular bridge to all-trans config- 
uration. Avalanche propagation is sustained by subsequent isomeric 
transitions of neighbouring molecules 5 . Negative avalanche pulses 
correspond to events of charging the surface traps with photo- excited 
holes, while positive ones correspond to their recombination with 
electrons in NPs, after the back transfer through the SAM linker 5 . 

Figure 2 shows examples of distributions of avalanches size, dura- 
tion, and average size vs. duration. The distributions follow power 
laws, suggesting that our system exhibits critical crackling dynamics. 
This is further confirmed by verification that experimentally deter- 
mined exponents fulfil the critical exponent relation 16 a — 1 = (t — 1) 
/ avz and, most significantly, that individual avalanches of various 
durations can be collapsed to a universal shape, using the scaling 
relation I(T,t) = T ll(JVZ ~ l f shape (tlT), as shown in Fig. 3. The collapse 
was obtained with exactly the same avz exponent as that of the 
dependence of the average avalanche size on duration. 

Avalanches in the investigated system clearly propagate in two 
dimensions. The values of obtained critical exponents are in extre- 
mely good agreement with the exponents obtained by numerical 
simulations of 2D random field Ising model 6 . The experimental 
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Figure 3 | Pulse shapes for BDT device before (a) and after (b) scaling 
collapse. Inset shows the skewness of the pulse shape as a function of the 
pulse duration. 
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exponent T exp = 1.48 ± 0.12, while the theoretical T th = 1.5. The 
exponent 1/gvz determined from the average avalanche size depend- 
ence on duration is l/(jvz exp = 1.76 ± 0.26, the theoretical = 
1.56, while the one calculated from the exponents relation is l/<JVZ ca \ c 
= 1.67 ± 0.29. Our results indicate that the molecular system may 
well belong to the universality class of the random field Ising model 
(RFIM) 1 " 3 . 

The pulse shape predicted by RFIM is a symmetric one. Our 
experimental collapsed pulse shapes, similar to Barkhausen noise 
experiments, are clearly asymmetric 1718 . Pulse asymmetry can be 
quantified through the pulse skewness y s = (i 3 /a 3 , where a is the 
standard deviation and \i 3 the third moment around the mean. As 
shown in the inset to Fig. 3, our pulse asymmetry is negligible for long 
and short pulses and peaks for medium length pulses around 350 ms. 

Since charges are less likely to tunnel through molecules in the 
gauche state, the disorder in 2D molecular layer can be associated 
with the amount of gauche defects in the SAM. In a very clean system 
the amount of such defects is low. Any avalanche will easily find a 
path of molecules in all-trans state across the entire system size. At 
the critical disorder, the avalanches will take all possible sizes and the 
distribution will fit the power low at all sizes. With disorder increas- 
ing above the critical point the percolation of avalanches will be 
hampered by increasing density of gauche defects, capable of halting 
the avalanche propagation, and causing the appearance of exponen- 
tial cutoffs in the power law distribution that eliminate large 
avalanches. 

We have observed that the size and frequency of occurrence of the 
avalanches decrease with illumination time, to the extent that after 
some 50-80 minutes avalanches practically disappear. Only by 
allowing the sample to rest in dark conditions for several hours, 
the avalanche patterns are fully restored. We ascribe time evolution 
of avalanches to the loading of FET surface states by holes transferred 
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Figure 4 | Time evolution of the avalanche size histograms, 
logarithmically binned and normalized, see the text, at different light 
intensities; high light intensity (a), intermediate light intensity using 
0.2 dB attenuator (b), and low light intensity using a 2 dB attenuator (c). 
Time records were taken at intervals of 16 minutes (0-16 min, 16-32 min, 
32-48 min, 48-64 min, 64-80 min). Arrows indicate the direction of 
elapsed time. 



from the photo-excited NCs. The probability of avalanche occur- 
rences decreases with increasing fraction of charge loaded surface 
states. When the amount N c of charges transferred to the surface is 
much higher than the amount N r of available surface traps, as in the 
case of high intensity light, the charge build-up is governed by surface 
properties through competition between charge capturing and 
releasing events. The accumulated charge Q(t) increases as Q(t) ~ 
N r {l — exp( — f/T e ff)}, where r e ff is the effective rate of charge build 
up, 1/Teff = 1/t c + l/r r , t c is the average rate of charge capture and T r 
the rate of charge release. For very low levels of the illumination 
power N c <C N r and charge build-up becomes controlled by the 
number of transferred charges N c , Q(t) ~ N C {1 — exp( — f/T eff )}, 
where N c (t) = N 0 {1 — exp(— yr)}, y depends on the tunnel barrier, 
and N 0 is the number of available photo-excited holes in NCs, what 
makes charging faster for higher level illumination. 

Alkene molecules in 2D SAM are coupled through Van der Waals 
forces associated with the intrinsic molecular electric dipole moment. 
The vector of the intrinsic moment of dithiol alkens self- assembled 
on the surface of GaAs points in the direction out of the surface 19 . 
Therefore, accumulation of photo-excited holes in the surface traps 
effectively decreases the moment and weakens the coupling. The 
change in the moment is proportional to the number of charges 
trapped on the surface. Since RJJ for a given system remains con- 
stant, weakened coupling decreases the critical disorder R c . Thus, at a 
given level of disorder in the system R > R c , the distance from the 
critical disorder will increase exponentially with time as (R — R c ) ~ 1 
— exp( — kt). 

Time evolution of avalanche size distributions at different light 
power levels of NDT molecules linked device is illustrated in Fig. 4. 
The arrows in the graph indicate the direction of elapsing time. In 
order to reduce the data scatter in the tails of the distributions we 
have binned the data using logarithmic binning. Moreover, for pur- 
poses of better comparison the distributions have been normalized in 

such a way that D(s)ds = 1. The major effect of elapsing illumina- 



tion time is the decrease in the size cut-offs through weakening the 
coupling / between the molecules. 

Theoretical models of disorder induced critical dynamics predict 
that the cut-offs S max in the avalanche size distributions should scale 
as the power law of the distance from the critical point S max ~ (R — 
Rc) ~ Va 13 - The time dependence of S max is illustrated in Fig. 5. Since R 
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Figure 5 | Time evolution of the size distribution cut-off S max as revealed 
from Fig. 4 for 0.2 dB and 0 dB attenuation of the light power. The time 
axis scale is sequent number of 16 min time record. Dashed lines are the 
best fits to the logarithmic time dependence of 5 max . The linear dependence 
implies that S max ~ (R — i? c )~ 1/CT . Inset: size distributions of avalanches in 
NDT and BDT linked devices under full power illumination. 
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— R c is proportional to {1 — exp( — /If)}, Fig. 5 actually presents S max 
as a function of ln^R — R c ). One may conclude that S max decays as 
the power law of the distance R — R c from the critical point. 
Unfortunately, at this stage we are unable to determine the critical 
exponent a from the data because we cannot determine the value of 
the parameter X. 

Lastly, the coupling between SAM molecules depends also on the 
molecule length, as testified by decreasing melting temperature of 2D 
molecular layers with decreasing molecular length 20 . Our experi- 
ments show that indeed S max for 1,9 nonanedithiol linker is much 
larger than S max for shorter and less strongly coupled 1,4 butane- 
dithiol linker, see the inset to Fig. 5. 

Discussion 

The growing interest in nonequilibrium phase transitions and 
advanced development of theoretical models of disorder-induced 
critical phenomena creates a need for experiments probing crackling 
noise dynamics in systems with changing disorder. If the model 
predictions for the connection between the avalanche statistics and 
the disorder can be proven in experiments, then future noise mea- 
surements can be used to extract information about the amount of 
disorder in a system, without the need for direct measurements of the 
disorder. These results amount to a new type of nondestructive test- 
ing. Most importantly, the predictions for the effects of disorder on 
the noise statistics are universal, i.e. they are expected to apply to 
many different systems. This implies that experiments on one of 
these systems can be used to predict the behavior in many other 
systems that belong to the same "universality class". 

The major body of experimental work concerning disorder 
induced criticality and associated crackling noise dynamics has been 
performed using magnetic systems exhibiting Barkhausen noise 2 . 
However, experimental verifications of effects of changing disorder 
on dynamic properties of crackling avalanches are still missing. In 
our paper we discuss dynamics of a hybrid organic-inorganic nano 
system that exhibits critical avalanche dynamics and show that the 
critical disorder can be experimentally tuned by light illumination. 
Our hybrid molecular system enables experimental studies of critical 
dynamics in which the distance from the critical point depends on 
density and length of linking molecules and can be tuned during the 
experiments by changing light intensity and illumination time. The 
system provides thus the tool box for in-situ studies of critical 
dynamics under changing disorder. The difference between the 
coupling of the BDT and NDT molecules leading to different levels 
of critical disorder in SAMs composed of such molecules present 
another means of controlled coupling. 

Therefore, the investigated hybrid system provides a flexible and 
controllable tool box for experimental studies of critical dynamics in 
which the distance from the critical point can be in-situ changed 
during the experiments. This dynamics is revealed by manifestations 
of critical charge transfer dynamics, observed through the conduc- 
tivity noise in the channel of field effect transistor to the surface of 
which the molecular layer links semiconducting nanocrystals illumi- 
nated by infrared light. The critical disorder in the molecular layer 
can be experimentally tuned by the duration and intensity of light 
illumination. In addition, the critical disorder can be adjusted by the 
length and density of molecules in the self-assembled layer, and tem- 
perature. The system can be easily fabricated using self-assembly 
methods. Lastly, our experiments show that charge transfer through 
alkyl monolayers, frequently used as elements of molecular electronic 
devices, is a complicated process that can evolve into unwanted noisy 
avalanche dynamics decimating the performance of practical devices. 

Methods 

The FET device. The FET was made by standard lithography techniques from a 
conducting n-GaAs layer, grown on a super-lattice buffer layer of 6 periods (Al, 
Ga)As, each period about 100 nm thick, and 5 nm thick insulating layer of GaAs 
covering the conducting layer. To increase the sensitivity of the device the conducting 



channel of the FET was shallow and localized between 20 and 50 nm from the surface. 
SAM linker of 1,9 nonanedithiol HS-(CH 2 ) 9 -SH (NDT), or 1,4 butanedithiol HS- 
(CH 2 ) 4 -SH (BDT) molecules is self-deposited between the source and drain FET 
electrodes from an absolute ethanol solution. When the device is illuminated by light, 
with energies smaller than the GaAs energy band gap and larger than the InAs energy 
gap, we excite only the NCs, leaving the GaAs parts of the transistor unaffected. When 
the NCs absorb light, the photo-excited holes are transferred to the FET surface states 
via the organic molecules, change the surface potential, and thus modify the 
conductivity of the 2DEG transistor channel. It should be emphasized that FET in our 
experimental arrangement operates not only as a light detector but also as an 
amplifier of the light effect. 

Nano crystals. We have used in this work InAs 6.5 nm core and InAs/CdSe/ZnSe 
core/shell/shell NCs with the same core size. The InAs nanocrystals were synthesized 
in trioctyl phosphine using InCl3 and TMS3 as precursors and underwent 
subsequently a size selection process. For the shell growth an ionic layer adsorption 
and reaction process has been used. The shell consisted of 1 monolayer (ML) of CdSe, 
followed by 4 ML of ZnSe. As a result, the emission from the NCs at 0.96 eV has 
improved by more than 2 orders of magnitude. 

Measurements. All measurements were performed in a specially designed rod and a 
cooler sitting inside a Faraday cage. We used a FEMTO DLPCA-200 Variable Gain 
Low Noise Current Amplifier that was located as close as possible to the rod contacts. 
The system was situated on a grounded optical table and the rod and amplifier were 
shielded and grounded to the same ground. Dual channel Stanford research SR785 
spectrum analyzer was used to measure the spectra and time response of the system. 
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